Hole-pair symmetry and excitations in the strong-coupling extended t — model: 

competition between d-wave and p-wave 
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We analytically calculate the ground state pairing symmetry and excitation spectra of two holes 
doped into the half-filled t — t' — t" — Jz model in the strong-coupling limit (Jz >> \t'\, \t"\). 
In leading order, this reduces to the t' — t" — Jz model, where there are regions of d-wave, s-wave, 
and (degenerate) p-wave symmetry. We find that the t — Jz model maps in lowest order onto the 
t' — t" — Jz model on the boundary between d and p symmetry, with a flat lower band in the 
pair excitation spectrum. In higher order, d-wave symmetry is selected from the lower pair band. 
However, we observe that the addition of the appropriate i' < and/or t" > 0, the signs of t' and 
t" found in the hole-doped cuprates, could drive the hole-pair symmetry to p-wave, implying the 
possibility of competition between p-wave and d-wave pair ground states. (An added t' > and/or 
t" < generally tend to promote d-wave symmetry.) We perturbatively construct an extended 
quasi-pair for the t — Jz model. In leading order, there are contributions from sites at a distance 
of \/2 lattice spacings apart; however, contributions from sites 2 lattice spacings apart, also of the 
same order, vanish identically. Finally, we compare our approach with analytic calculations for a 
2x2 plaquette and with existing numerical work, and discuss possible relevance to the physical 
parameter regime. 
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PACS numbers: 71.27-|-a,74.20.Mn 

I. INTRODUCTION 

In recent years a number of experiments, particularly 
phase-sensitive ones, have indicated that the pair symme- 
try of hole-doped cuprate superconductors is at least pre- 
dominantly cia;2_j^aiiSii. Theoretical and numerical stud- 
ies of the two-dimensional Hubbard, t — J, and related 
models have also suggested d^2_yi pairing^iSiSiSiSiS, and 
studies of Hubbard and t — J models on 2 x 2 plaquettes 
have provided an intuitive picture of how d-wave sym- 
metry might ariseiSiii. However, there are few rigorous 
theoretical results in this general area. 

Different experimental techniques have indicated that 
a pseudogap with the same symmetry as the su- 
perconducting gap persists above Tc in undcrdoped 
cuprates^- ^^•^ ^d'^d^d? . This, along with the short high- 
Tc coherence length^, is qualitatively consistent with 
a strong-coupling picture, where pairs can preform at 
T > Tciii. Numerical work has in addition suggested 
that the t — J and t — models have many similar 
propertie o'^d^d^i^° (see however Ref. and that the 

t — Jz model may hence provide a suitable starting point 
for understanding t — J behavior—. 

Reductions from Cu02 three-band and similar 
modelsS^, as well as comparison with ARPES results for 
a single doped holeS4, suggest that, besides a nearest- 
neighbor (NN) i, the next-NN (NNN) t' and next-NNN 
(NNNN) t" hoppings may also be substantial. Nu- 
merical calculations and theory have explored some of 
the qualitative effect that the sign of t' has upon hole- 
pairingS^iSSiSi. Given these above results, it is interesting 
to try to gain a better understanding of the pairing prop- 



erties of the extended t— J model. In addition, performing 
such a study analytically could help provide guidance to 
future numerical work aimed at exploring specific regions 
of the model's phase diagram. 

To that end, we consider in this paper two holes doped 
into the half- filled t — t' — t" — Jz model in the strong cou- 
pling limit {Jz >> \t\, \t'\, \t"\). We calculate the symme- 
try of the hole pair in the ground state as well as the pair 
excitation spectrum. We do not explicitly consider the 
issues of phase separation or whether superconductivity 
actually occurs. We consider first the t' — Jz model, and 
show how singlet pairs can be constructed from our solu- 
tions. We next discuss the t' — t" — Jz model, and then 
the t— Jz and t — t' — t" — Jz models. For the t—Jz model, 
we perturbatively construct an extended quasi-pair. As 
a step towards exploring the range of validity of our ap- 
proach, we compare with results for a 2x2 plaquette and 
with numerical studies. Lastly, we discuss implications of 
our results for the physically relevant parameter regime, 
among which a phase competition scenario is one possi- 
bility. 

Specifically, we consider the Hamiltonian 



H — Hq + Hi + H2 + H3 



(1) 



where 



Ho — Jz X^K^y {('^a:,i;'^a;-|-l,i; + ^x,y^x,y+l) 



T^^x.y^x + l.y ~^ '^x.y'^x,y-\-l^} ^ (2) 



Hi = (-t) E 



^ - J- H.c.) 

+ {cl,y,a^x,y + l.a + H.C.)}, (3) 



x,y,a {i^x,y,a'^-c+i,y.<^ 



2 



H. 



2 = i-t') J2x.,y.a- {(4,y,<TC2: + l,y+l,£T + ff-' 



+ H.C.)}, (4) 



and 



+ {cl.y,aCx,y+2,a+H.C.)}. (5) 

Here, x and y denote the coordinates of an LxL lattice 
with periodic boundary conditions and even L (with, ex- 
cept for the t' — Jz model, L > 2), and a = ±1 (t,i) 
refers to electron spin. Cx^y.a = Cx.y.ai^ ~ "-x.j/.-cr)! en- 
forcing the condition of no double occupancy. S^ y = 
1/2 {nx,y,-\ - nx,y,i) and Ux^y = n^.y,] + "x.^a- We do 
not explicitly consider here the spin-flip part of the mag- 
netic interaction 



(t)i: 

x,y 



{{^x.y'^x+l.y + ^x,ySx,y+l) + H.C.^ ,(6) 



where 5'+,, = ct 



x,y 



Zx,y,i and S^^y = cl^y^^Cx,y.,]. (The 
full t — t' — t" — J model is recovered when — Jz ) 

At half filling each site is occupied by exactly one 
electron, and the doubly degenerate ground state of 
Hq is then that of a Neel antiferromagnet. We choose 
|$Q > to denote the state with electron spins a{x,y) = 
(— 1)^+^ and |$h > to denote the state with a(x,y) = 
(— 1)^+*'+^. We define the operator a^^y — Cx^y^a{x,y) with 
(7{x,y) = (-1)^+'', and the operator bx,y = c^,y,o{x,y) 
with a{x,y) = (— 1)^+^+1 Although our calculations and 
results are independent of the ordering convention cho- 
sen, we will denote for specificity 

>= (a^...al^i)...(a[^2-4,2)(aI,i-«M)|0 >' (7) 

where |0) is the state with no electrons, with an analogous 
definition for |$b >. 

We now dope the half- filled state |$a > with two 
holes and consider the strong-coupling limit (Jz >> 
\t\,\t'\,\t"\). In this limit, there will be an energy cost 
of order Jz if the two holes are not NN. Hence, to zeroth 
order, the (highly degenerate) two-hole ground state is 
spanned by the set of all NN hole pairs. We denote the 
state with a horizontal NN hole pair at sites [x, y) and 
{x + \,y) as 



\^x,y ^ — ^: 



x+l,y^x,y 



ax,v\'^a >, 



(8) 



and the state with a vertical NN hole pair at sites (x, y) 
and [x, y + 1) 



\Vx,y >= ax,y+iax,y\^a > 



(9) 



The \hx.y >'s and iv. 



x,y >'s provide a complete, orthonor- 
mal basis for the two-hole ground state of Hq correspond- 
ing to |$a >■ 

It costs an energy of order Jz if one of the NN holes 
hops to a NN site through the hybridization matrix ele- 
ment t. However, there is no energy cost for hops corre- 
sponding to t' or t" , as long as the two holes remain NN 



after the hop. Thus, to lowest order in 1/Jz, it is only 
necessary to diagonalize the Hamiltonian H2 + H3 in the 
subspace spanned by the \hx.y >'s and \vx,y >'s; i.e., it is 
only necessary to consider the t' — t" — Jz model. We note 
that in this limit the t' — t" — Jz model becomes isomor- 
phic to the strong-coup ling limit of the antiferromagnetic 
van Hove model of ref. 123. 



II. t' MODEL 

We consider first the t' — Jz model, involving only the 
H2 (diagonal) hopping term. Defining 



and 



^ e 



(10) 



(11) 



with kx,ky = 0, 1...L— 1, we obtain the lowest order 
wave functions 



> 



with energies 



± 



V2 
± sgn(t') e- 



\vk^,k,>] (12) 



(13) 



Since < sm[T:kx/ L),sm{'Kky/ L) < 1, the minus sign 
gives the branch of lower energy. The lowest energy state 
I'f/'o"'' >j with energy — 4|t'|, occurs when kx = L/2 and 
ky — L/2 (i.e., (tt, tt)). Rewriting in terms of the ax,y's 
and neglecting overall phase factors, one obtains 



1 



lV2 



T,x,y {-^T'^^{ax+Uyax,y 



sgn(i') }|<i>a >(14) 



When t' > (sgn(i') = 1), the sum over hole pair oper- 
ators in Eq. 1141 changes sign upon a 90 degree rotation 
around a lattice point, giving the pair d-wave symme- 
try (specifically, dx2_y2mi^). When t' < 0, there are no 
such sign changes, giving s-wave symmetry (specifically, 
extended- s2ii). 

If one adds to Eq. the appropriately-phased pair 
operator for two holes doped into the ground state >, 
given in 

>= ^ Ex,y i-ir+'+Hbx+l.ybx.y 

~-Sgn{t')b^,y + lb^^y}\<i>b >,(15) 



3 



one obtains for t' > the usual NN singlet dx'2_y2 pair 
operator 



2J^ '^x.y \{.'^x,y,]Cx+l,y,l Cx,y,lCx+l,y,]) 



-{cx,y,-\Cx,y+l,i - Cx,y,iCx.,y + l,{)^ 

with t' < giving the analogous singlet extended-s oper- 
ator 

A3 = 

^ Tlix,y \^{^-x,yAC-x + l,y,i — Cx^y^lCx + l,y,-\) 

+ (cx,y,TCa;,y+l,i " Cx ,y ,iCx,y+lA)'^ ■i^^) 

With different relative phases, one can also obtain types 
of d-wave or s-wave to = triplet pairs; because quantum 
spin fluctuations are not included in the t— Jz model, the 
cases cannot be differentiated at this level. 

One can better understand the dependence of the 
t' — Jz (and t' — J) pair symmetry on sgn(f') by con- 
sidering phase transformations of electron creation and 
destruction operators. Specifically, we consider transfor- 
mations of the form Cj — e^^^^^dj, where j is some gener- 
alized coordinate referring to both orbital and spin and 
9{j) is some function of j. 

First, as background, let P denote some arbitrary 
product of electron creation and destruction operators, 
some of which operators may be the same, with coordi- 
nates referring to orthogonal states. For example, for a 
one-dimensional chain with one orbital per site, one could 
have P = c\^c\^ciicl^c\y 

As previously, let |0) denote a state with no electrons. 
Then, (0|P^P|0) if it is nonzero reduces to 

(0|ptp|0) = (0\l[c,c]\0) = (0|n(l-^^.)|0> (18) 



for some subset {j} of the generalized j coordinates. 



Hence, any phase transformation of the form cj 



does not affect the value of (0|PTp|o). 

Now, consider (0|P/Pi'|0), where Pi and P;/ are two 
different products of creation and destruction operators. 
Each nonzero (0|P/P;' |0) also reduces to the general form 
(0| rijj} Cjcj |0), which is again unaffected by the trans- 
formations Cj = e^^^^'^dj. Lastly, consider the expecta- 
tion value of an operator O = a;P; in the (possibly 
unnormalized) state biPi\0), where the ai and bi are 
coefficients. Then, 



(O) 



ai"blbi,{0\P^PrP,\0) 
Eu'btbi,{0\PlPi,\0) 



(19) 



Again, each term of the numerator and denominator 
reduces to the form (0| Hfj} ^jc] |0), which is invari- 
ant under cj — e^^'^^^dj. Hence, an arbitrary phase 
transformation on electron creation and destruction op- 
erators has no effect on fermion operator expectation 
values in fermion states. Also, since an electron cre- 
ation/destruction operator referring to a particular ba- 
sis can always be written as a sum of operators referring 
to a different orthogonal basis, it is not even necessary 
that the creation/destruction operators under considera- 
tion refer to a particular orthogonal basis, as long as the 
phase transformations are consistent. 

Now, we wish to find a particular phase transformation 
which is equivalent to reversing the sign of t' . Temporar- 
ily dropping the spin coordinate, we choose for simplicity 
a transformation of the form 



„ _ Jiip^x+ipyy) • 
^x^y — iLx 



(20) 



where (jjx and 4>y are constants. (We do not include a 
constant phase factor as it has no relevant effect.) We 
require: 

cl+^^y+iCx,y + H.c. = {-l)[dl^^y^^dx,y + H.C.] (21) 
and 

cl_^.y + lCx,y + H.C.= {-l)[dl_^y^^dx,y+H.C.]. (22) 

Solving the previous equations, one can obtain 

0, -I- (j)y = {2p + l)Tr (23) 

and 

- (j>y = {2q + l)7r (24) 
for arbitrary integers p and of which solutions are 

Ck.j/ = ( 1) d,x^y, (25) 

or 



i-l^dx 



(26) 



Restoring the spin coordinates leads to the four following 
transformations equivalent to changing the sign of t': 



^x,y, 



i-lTdx,y,h Cx,y,i = {-\Ydx,y,i\ (27) 



Ca;,y,T = (^ 1 T ; Cx,y,i = i^^)^ dx ,y ,1', (28) 

'^x,y,1 — ( 1) dx,yA, ^x,y,l — ( dx^y^, (29) 

Cx,yA = {~iydx,y,^, Cx.y,i = {-Vfdx^y^. (30) 

None of the above four transformations changes the 
sign of the Jz term (Eq. However, the transforma- 
tions of Eqs. and 1201 (as well as more general trans- 
formations not discussed here which are different on dif- 
ferent sublattices) do not leave the i/^ of Eq. El and 



hence the spin-spin interaction of the t' — J model, in- 
variant. The transformations of Eqs. 1271 and which 
do leave H±^ invariant, change the dx2_y2 symmetry of 
the pairs of Eq. 1141 to s and the d^2_yi singlet pairs 
of Eq. El to extended-s singlets and vice- versa. Hence, 
dx2^y2 singlets of the form of Eq. 1161 are transformed to 
extended-s singlets when the sign of t' is changed in the 
t' — J model, and vice-versa. More generally, operators 
are transformed according to Eq. 1271 or E51 when the sign 
of t' is changed. (The two transformations give the same 
results for operators with the appropriate symmetries, 
which generally includes the operators of interest). This 
is different from ref. |2^ However, the main arguments 
and conclusions of that paper remain unchanged. 



t' 



III. 



t' - 1" 



J, MODEL 



For the more general t' — t" — Jz model, one obtains in 
lowest order the (unnormalized) wave functions 



wikyy 



(2t"){s: 



^ -S^)±T 



with energies 



+ 



where Sx = sin(7rfca;/L), Sy = sn\{nky / L) , and 



(32) 



(33) 



As a function of t' and t", we find that the ground state 
symmetry of the pair is as shown in Fig. 1. 

The s-wave and d-wave operators are of the form in 
Eq. El The ground state associated with p-wave pairs 
is, in contrast, highly degenerate. The multiple p-wave 
pair operators can be either px 

1 ^ — ^ 2-Kikyy 

\^ (ky) = J—^ 2^ e ^ax,y{ax+uy - ax-uy) (34) 



x,y 



or py 



lV2 



0'x^y{cix,y+l ^a:,y— l)- (35) 



In leading order, the Px states have energies independent 
of ky, and the Py states have energies independent of 
kx- Both p-wave pair operators change sign under a 180 
degree rotation. 



IV. t - J, MODEL 

We next consider the strong-coupling limit of the t—J^ 
model. To lowest order, we find that this maps onto the 



FIG. 1: Hole pair symmetry in the strong-coupling limit of 
the t' — t" — Jz model as a function of t' and t" . The legend 
is as follows: "D" denotes d^2_y2, "P" denotes px or py, and 
"5"' denotes extended-s. The axes cross ai t' = and t" — 0. 



above strong-coupling limit of the t' — t" — Jz model with 

2 / t2 ■ 

3 IX 



h,eff. - tt.eff. 



(36) 



From Eq. 1321 the lower band of the pair excitation spec- 
trum then becomes flat, with wave functions 



Sy\hk^,ky > 



e ^ Sx\Vk^ky > 



}• (37) 



Flat pair bands were also found^^'^" for related (though 
different) models and/or treatments. In ref. [^J, a five- 
fold degeneracy of strong-coupling t — Jz pairs of pure d 
or p symmetry was noted. 

We see from Eq. (SHI that, to lowest order, the strong- 
coupling t — Jz model lies on the (rightmost) boundary in 
Fig. 1 between d-wave and p-wave symmetry, providing a 
simple picture for competition between these two states. 
In the next higher order, neglecting constant additive 
terms, the energies of the lower pair band separate into 



{cl +cl+ ACxCy - ilCx 



31cy + 56},(38) 



where here Cx — cos(27rfca;/L) and Cy — cos(27rfcj,/L). 
Two-hole lower band dispersion curves were previously 
calculated numerically using a variational method^^ and 
series expansionsSS. We then find (in agreement with 
refs. I23 and |2^ that the pure d-wave {t' > 0) state of 
Eq. ^^is selected as the lower pair band ground state. 
However, the closeness to p-wave symmetry may provide 
an explanation for the low-energy p-wave "quasi-pair" 
peaks seen numerically in small t—J and t—Jz clustersi^. 

We also note that several different techniques have sug- 
gested that the symmetry of a doped hole pair in the 



5 



t— Jz and/or t — J models may be p-wave for some range 
of generally intermediate or small J^/t or j/^20,3i^32,33_ 
However, this has not been rigorously confirmed one way 
or the otheri^i^. 

One can perturbatively construct increasingly ex- 
tended quasi-pair states for the t — model. Combining 
results for the NN d-wave pair operators for ground states 
|$a > and |(f>f, >, one finds the lowest order correction 
for the singlet pair operator of Eq. 1161 

--H)(x)(i)S' 

x,y,(T 

I {^x+l,y,a'^x+l,y,-a^ Cx,y+l,-(T'^x,y+l,iT) 
+ 1 .y + 1 .a^x ^y ^(T 
+ {'^x+l,y,<y'^-x + l,y,-a^ '^\:.y-l,-a^x,y-l,cr) 

^x-\-l,y — l,a^x,y,a^- ('^^) 

When operating on the appropriate Neel state, each of 
the above terms consists of a diagonal hole pair "dressed" 
with a singlet pair of electrons straddling the bond con- 
necting the pair of holes, as was found in numerical t — J 
simulations^"'^. This form basically arises from the disrup- 
tion in the local spin order when an electron of particular 
spin orientation and one hole of a NN hole pair exchange 
sites to form a diagonal hole pair. We note that the con- 
tribution from pairs a distance of two lattice sites apart, 
nominally also of order t/Jz, vanishes identically in this 
order. (Variational calculations found a reduction of such 
a contribution"^^.) This vanishing may provide an expla- 
nation for why only NN and diagonal hole correlations 
appear to dominate in the t — J model near half filling 
for moderate to large J/-^L^. 

If one adds the necessary terms to the operator of Eq. 
1391 to impose a type of spin rotational invariance, one 
obtains the composite pair operator invented in ref. [s^ 
to give a diagonal singlet pair with (1^2 _y2 symmetry. 
The type of operator of ref. QJi was also generatedli in a 
2x2 plaquette by applying e~'^^ to the pair operator of 
Eq. ^1 The operator of Eq. 123 which emerges naturally 
from perturbation theory, also has dx2_y2 symmetry. 

We also note that, since we calculate energy spectra 
and wave functions, our results and approach can be used 
to calculate finite-temperature and real frequency prop- 
erties. Our results can also be easily extended to periodic 
n-leg ladders with even n. However, we do not pursue 
those issues here. 

We in addition note that, in contrast to Eqs. [771 
through 1301 the phase transformation that is equivalent 
to changing the sign of t is given by 



(40) 



This transformation leaves the spin-spin interactions Hq 
and the H± of Eq. El invariant. Under this transforma- 
tion the symmetries of the hole pair operators of Eqs. 1141 
1161 and 1391 all remain the same. 



V. t-t' - t" - J, MODEL AND INTERMEDIATE 
COUPLING 

For \tl \t'\, \t"\ « Jz; \t\ » \t'\,\t"\; and \t'/Jz\ > 
{t/ JzY or \t" / Jz\ 3> {t/ JzY (so that lowest order terms 
dominate), one can use Eq. 1361 to define a 



eff. 



eff. 



and 



eff. 



, 2 e 

= t H 

3 J. 



„ 2 e 



(41) 



(42) 



One can then simply use Fig. 1 as a guide to pair sym- 
metry. 

As one step towards investigating the further range of 
validity of our approach, we performed analytic calcula- 
tions of the t — t' — Jz — J± model (sec Eq. (HI on a 2 x 2 
plaquette. In what follows, the four plaquette sites are 
consecutively numbered as one goes around the plaquette 
edge. Spin-spin interactions and t hoppings are between 
consecutive sites, t' hoppings are between diagonal sites, 
and there are no periodic boundary conditions (though 
such boundary conditions would only renormalize param- 
eter values). 

With no holes (four electrons, one per site), we found 
that the ground state had dx2-y2 symmetry, as is the case 
for the t — J model on a plaquett»i2*ii, for all Jz > 0, 
J_\_ > 0. (There are two degenerate states only when 
J± = 0, both of d-wave symmetry). 

For two holes, when 

1 



t' > <^,„,,. = - [{Jz + Ji) - vIjT+jIW+g^] (43) 

(the more "physical" parameter regime), the ground 
state energy is 



1, 



where 



(44) 



(45) 



(We assume Jz > 0, J± > 0, and cither t or t' ^ 0). 
Dcfiningii 



(46) 



with Aj^ = Aj-, where i and j refer to plaquette sites. 



t. = At. 

the ground state is a linear combination of the two states 
1 



|^i) = -^(Al3 + Ay|0) 



(47) 



and 



1^2) = :.(AL + aL + Al + Al)\0). (48) 
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The ground state hence has s-wave symmetry, as was 
earher found for the two-hole doped plaquette t — J 
modelifliil. As discussed prevfousljiiSiii, this imphes 
that the hole pair symmetry is d-wave. We note that 
with i' = 0, which automatically satisfies Eq. 021 the 
hole pair symmetry for the plaquette t ^ Jz — Ji_ model 
is always d— wave, as in the strong-coupling limit of the 
general t — model. (However, because there are no 
lattice sites a distance of two lattice spacings apart on 
the plaquette, for t/J^ « 1 the plaquette t—J^ model 
reduces to a i' — model with t' > rather than a full 
t' — t" — Jz model, enhancing d-wave pairing compared 
to a larger lattice.) 

However, there is a level crossing at t' = t'cross.^ ^-i^d 
when t' < t'^ross. the ground state has energy 

Eo = ^^-^ + 2t' (49) 

and wave function 

|V'3) = ^(AL-At3 + Aj4-Ali)|0). (50) 

IV'a) has d-wave symmetry, implying that the hole pair 
symmetry is now s-wave. We note that with t = 0, giving 
^'cross — the t' — Jz — J± model hole pair symmetry 
on a plaquette is d-wave for t' > and s-wave for i' < 0, 
as in the general t' — Jz model strong coupling case. 

As mentioned previously, some approaches have sug- 
gested p-wave hole-pair symmetry for some range of gen- 
erally intermediate or small J/t or Jzf&^ii^^i^. How- 
ever, other t— J numerical results^Si^^iSiSi indicate that 
d-wave features of the strong-coupling limit may persist 
down to physically relevant intermediate coupling {J/t « 
0.3 — 0.5). This and the above plaquette results are con- 
sistent with the strong-coupling t — t' — t" — Jz model 
as a potentially useful starting point to explore the more 
physically relevant intermediate-coupling t — t' — t" — J 
model. 

As a rough guide to the general effects of t' and t" on 
pairing symmetry, one can begin with a d-wave state in 
Fig. 1 arising from the t term, presumably in the upper 
right quadrant. Assuming pairs of pure symmetry and 
that strong coupHng qualitatively extends to \t'\, \t"\ ~ 
J, Jz, both t' > and t" < will tend to move one deeper 
into the d-wave region; however, t' < and t" > will 
tend to move one towards, and perhaps into, the p-wave 
region. Previous t — J work has argued and/or indicated 
numerically that t' > strengthens d-wave pairing while 
t' < weakens i1i2Si21, and it was found on a 32-site lattice 
that a particular t' < and t" > together favored a p- 
wave pair—. Also, if one starts with a possible t ~ Jz 
ov t — J p-wave pair (upper right quadrant), i' < and 
t" > would tend to move one deeper into the p-wave 
region while i' > and t" < would tend to move one 
toward, and perhaps even into, the d-wave region. 

Based on this and our strong-coupling results (assum- 
ing pairs of pure symmetry), we show in Fig. 2 qualita- 
tive predictions of the hole pair symmetry for the t—t'—Jz 
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FIG. 2: Qualitative diagram of proposed hole pair symmetry 
for the t — t' — Jz model. "D" , "P" , and "S"' denote the same 
as in Fig. 1, and no prediction is made for region "[/". Axes 
cross at t' /t — and intermediate J^/t (~ 0.3 — 0.5). The 
possibility of p-wave symmetry for t' — and intermediate 
Jz/t is not shown. 



model. We do not show in Fig. 2 the possibility of a 
cross-over to p-wave symmetry for the t — Jz model men- 
tioned above. We believe the predictions shown apply to 
the t — t' — J model as well, with a comparatively smaller 
p-wave region due to larger energy differences between 
t — J p-wave and d-wave pair statesi^. An additional 
t" > would tend to enlarge the p-wave region, and an 
additional t" < would tend to enlarge the d-wave re- 
gion. Note that in Fig. 2, the horizontal axis cuts the 
vertical axis at intermediate Jz/t (~ 0.3 — 0.5). 

t > for the hole-doped cuprates, and estimates for t' 
and t" are typically in the ranges t' « (— 0.1)t — {—0.5)t, 
and t" « 0.0- (0.3)t, generally with \t'\ > \t"\ Both 
these signs of t' and t" are those which could tend to drive 
the pair symmetry to p-wave, raising the issue of the hole 
pair symmetry in the intermediate-coupling regime. We 
note that s-wave is also possible, though we believe it less 
likely at intermediate coupling. 

It would be interesting to try to determine numeri- 
cally whether the symmetry of two doped holes in the 
t—t'—t" — J model is in fact d-wave for the experimentally 
relevant values of t,t',t", and J (e.g., J/t « 0.3 — 0.5). 
However, drawing conclusions from exact diagonalization 
and other current app roaches may be challenging (see, 
e.g., refs. |2^ and l33t) due to the possibility of uncon- 
trolled finite-size or other errors, and such approaches 
sometimes give confiicting results. If the symmetry were 
established to be p-wave rather than d-wave, it would 
suggest that the t — t' — t" — J model by itself could be 
incomplete as a model for high- Tc superconductivity. In 
that case, one possibility for restoring d-wave symmetry 
could be the addition of electron-phonon coupling in the 
d-channelS. In either case, it may also be of interest to 
explore how the existence of or nearness to p-wave sym- 
metry, which effectively reduces the dimensionality of the 
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hole pair wave function from 2D to ID, would correlate 
with possible stripe phases observed experimentallji22ii£ 
and in numerical simulation3^i*^Siii. 

Another interesting possibility is the existence of 
a novel exotic state which may be close in energy 
to the dx2-y2 superconducting cuprate ground state. 
Such a competing phases scenario'*^ could become even 
more involved if one takes in account recent numerical 
calculations^ pointing to the existence of strong ferro- 
magnetic fluctuations in the vicinity of doped holes for 
realistic values of i' < 0. 



VI. SUMMARY 

In summary, we have investigated analytically the 
ground state pair symmetry and excitation spectra of two 
holes doped into the half-filled t — t' — t" — model in 
the strong-coupling limit. In lowest order, this reduces to 
considering the t' — t" — Jz model, where we found regions 
of ground state d-wave, s-wave, and (degenerate) p-wave 
symmetry, depending upon the signs and relative magni- 
tudes of t' and t" . We next found that the t — model 
in lowest order was on the boundary between c?-wave and 
p-wave pair symmetry, with a flat lower pair dispersion, 
providing a simple picture for the competition between d 
and p symmetries. In higher order, d-wave symmetry was 
selected from the lower pair band. However, because of 
the closeness to p-wave symmetry, we predict that the ap- 
propriate t' < and/or t" > added to the t—Jz oi t—J 
models with intermediate to large Jz or J should drive 
them into p-wave pairs, and perhaps even p-wave super- 



conductivity. These signs of t' and t" are those found 
in the hole-doped cuprates. (In contrast, t' > and/or 
t" < tend to promote d-wave symmetry.) This p-wave 
tendency could be strengthened following results which 
suggest p-wave pair symmetry for intermediate or small 
J/t or Jz/t in the t - J or t - Jz modela^°i'^^i^^i^^, though 
such results have not been rigorously confirmed^^i^. 

We constructed a perturbative correction to the 
nearest-neighbor d-wave pair, giving a more extended 
quasi-pair, and found that it was similar to the d-wave 
composite operator invented in ref. [s^ and qualitatively 
consistent with previous numerical resultsii. The quasi- 
pair included a contribution from sites a/2 lattice spac- 
ings apart, but the same-order contribution from sites 2 
lattice spacings apart vanished identically. The structure 
of the quasi-pair derived from the disruption in local spin 
order under the exchange of one of the nearest neighbor 
holes and an electron. 

We explored ranges of validity of the perturbative ap- 
proach of this paper using a 2 x 2 plaquette and re- 
sults from other -^yorkSSi^^i^Si^i. Lastly, we discussed 
implications for the experimentally relevant parameter 
regime. These included the possibility of p- wave symme- 
try for two doped holes, which would suggest that the 
t — t' — t" ~ J model could be incomplete as a high- Tc 
model, or perhaps a phase competition scenario between 
d-wave superconductivity and a p-wave state. 
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